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We study the Kondo effect in muitipie-dot systems for which the inter- as weii as intra-dot 
Coulomb repulsions are strong, and the inter-dot tunneling is small. The application of the Ward- 
Takahashi identity to the inter-dot dynamical susceptibility enables us to analytically calculate 
the conductance for a double-dot system by using the Bethe-ansatz exact solution of the SU(4) 
impurity Anderson model. It is clarified how the inter-dot Kondo effect enhances or suppresses 
the conductance under the control of the gate voltage and the magnetic field. We then extend 
our analysis to multiple-dot systems including more than two dots, and discuss their characteristic 
transport properties by taking a triple-dot system as an example. 
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I. INTRODUCTION 

Recent advances in semiconductor processing have 
made it possible to fabricate various nanoscale materi- 
als with tunable quantum parameters, revealing various 
aspects of quantum mechanics. Quantum dot [lj, |2| is 
one of the interesting nanoscale materials. In particular, 
a lot of works on the Kondo effect in single quantum dot 
systems have been done both theoretically and experi- 
mentally [3, 0, la, lI 0, la 111 lii! EH LL=f • More recently, 
double-dot systems or s yste ms with more than two dots 
have been investigated JT3I IT3. Iil| . In this connection, 
the Kondo effect in double-dot systems have been studied 
intensively [H [H H |H E| . 

Most of the previous studies on multiple-dot sys- 
tems have treated the intra-dot Coulomb repulsion but 
have ignored the Coulomb repulsion between quantum 
dots (inter-dot Coulomb repulsion). We especially fo- 
cus on the effect of the inter-dot Coulomb repulsion here 
H3- E3- EH] and study how such electron correlations af- 
fect transport properties. Recently, Borda et al. have in- 
vestigated properties of the Kondo effect in such double- 
dot systems with a magnetic field by the numerical renor- 
malization group method |24j . which may explain the 
Kondo effect observed experimentally by Wilhclm et al. 
in a double-dot system 22] . A remarkable point in the 
above double-dot systems with inter-dot Coulomb repul- 
sion is that enhanced charge fluctuations between the 
quantum dots induce the "inter-dot Kondo effect" , which 
plays an important role to determine transport proper- 
ties of the systems. Since the inter-dot Kondo effect is 
caused not by spin fluctuations but by charge fluctua- 
tions between two dots, its influence appears significantly 
when the dots are connected in series. In particular, by 
changing the gate voltage or the magnetic field, one can 
control the conductance via the inter-dot Kondo effect. 

In this paper, we investigate transport properties of 
the double-dot systems with strong intra- and inter-dot 
Coulomb repulsions mentioned above. We exploit a novel 
method to treat the Kondo effect at absolute zero: the 



application of the Ward-Takahashi identity enables us to 
use the Bethe-ansatz exact solution of the SU(4) impurity 
Anderson model to our double-dot system. Our calcula- 
tion clearly shows that the inter-dot Kondo effect plays 
an important role on transport, which can be controlled 
by the gate voltage and the magnetic field. Our method 
is also applicable to multiple-dot systems including more 
than two dots. We explore the Kondo effect in such dot 
systems by taking a triple-dot system as an example. 

This paper is organized as follows. In section |H1 we in- 
troduce the model Hamiltonian and outline the method 
to treat our double-dot system: how the Bethe-ansatz 
exact solution can be used to compute the conductance 
at absolute zero. In section ITTT1 we discuss the results for 
the conductance with particular emphasis on the gate- 
voltage control and the magnetic-field control. In section 
IIVI we extend our method to a triple-dot system, and 
discuss its transport properties on the basis of the exact 
solution. We also mention how we can treat generalized 
multiple-dot systems including more dots. A brief sum- 
mary is given in section Ivl 



II. MODEL AND METHOD 

We describe our model and method by taking a double- 
dot system connected in series, which was proposed by 
Borda et al. |24|. The setup is schematically shown in 
FIG. n where not only the ordinary Coulomb repulsion 
U, which works inside each dot, but also U' between the 
dots are introduced. We assume that the inter-dot tun- 
neling t is small and the gate voltages are such that the 
lowest-lying charged states are restricted to the config- 
urations of singly-occupied states, (nn, ni,) — (1,0) and 
(0, 1), where is the number of extra electrons on 

the right (left) dot. This situation is realized in the con- 
dition, 



\E(1, 0)|, \E(0, 1)| <E/,f/' 



(1) 
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where E(nR,riL) is the energy level in the dots mea- 
sured from the common chemical potential of the two 
leads. The states (1,0) and (0, 1) have a spin S = 1/2, 
associated with the extra electron on the double dots. 
Then at energies below the charging energy of the dou- 
ble dots, dynamics of the double dots is restricted to 
the subspace with the 4 possible configurations of {S z = 
±1/2; tlji — til = ±1} in addition to the unoccupied state 
of n R = n L = 0. 

The above double-dot system, in which both of intra- 
and inter-dot Coulomb repulsions are sufficiently strong, 
may be modeled by the highly correlated degenerate An- 
derson Hamiltonian Ha {U,U' — > oo) supplemented by a 
inter-dot tunneling term Ht, 
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H A = /^4r(^|^w(z)+E £ * 

£T,T <J,T 

+ VJ2 f dxS(x)[\aT){0\c aT (x) + h.c], (2) 

(TT ^ 

H T =t K)(<tt'| (3) 

<t,t^t' 

where c] 7T (x)(c crT (x)) creates (annihilates) a conduction 
electron at a position x with spin a(= ±1/2) and "orbital 
index" r. Here we have represented conduction electrons 
in the leads in the low-energy continuum limit by assum- 
ing that its density of states is constant, 1/27T. Also we 
have introduced the orbital index t — 1/2 (-1/2) to spec- 
ify an electron occupying the left (right) lead, which is 
also used to label the left (right) dot. A state \ar) in the 
double dots located at x = denotes a singly occupied 
state and |0) denotes an unoccupied state. 

We will discuss transport properties of the system un- 
der the gate- voltage control or the magnetic-field control. 
It is thus convenient to write down each energy level e dT 
as, 



e a dT = e + 6E ■ t + E. 



Z ■ cr, 



(4) 



where SE (Ez) is the energy difference between the two 
dots (Zeeman energy). Note that the system possesses 
SU(4) symmetry with respect to spin and orbital degrees 
of freedom at 8E = Ez = 0. 

We note here that the Bethe-ansatz exact solution 
can be obtained for the above four-component Anderson 



Hamiltonian Ha [25|, which is referred to as the SU(4) 
Anderson model henceforth. However, this method al- 
lows us to calculate only static quantities, so that we 
cannot apply the exact solution to transport quantities 
straightforwardly. In the following, we outline how we 
can overcome this difficulty to calculate the conductance. 

Let us begin with the expression for the conductance in 
the above double-dot system connected in series, which 
is obtained in the second order in the tunneling Hamil- 
tonian Ht between two dots, 
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where \ops{^) is the analytic continuation (iio n — ► w+iO) 
of the dynamical "orbital pseudo-spin" susceptibility for 
the SU(4) Anderson Hamiltonian Ha (without Ht), 

fP . 
Xo P s(iu n )= dre*"" < TT + (r)T_(0) >, (6) 
Jo 

with the time-ordering operator T. The corresponding 
SU(2) operators are defined as 

1 



T ? = i^{nR-n L ), 



f ± = £|<7±1/2)<<7T1/2| 



(7) 



These orbital pseudo-spin operators properly describe 
inter-dot charge fluctuations. As defined above, the 
eigenvalue r = ±1/2 of T z specifies which dot an electron 
occupies. Eq. (JSJl means that the low- frequency inter-dot 
"orbital" susceptibility is essential to determine the con- 
ductance. 

Although the low-frequency susceptibility is difficult 
to calculate in general, we can make use of sophisticated 
techniques developed in the study of the NMR relax- 
ation rate in dilute magnetic alloys |26l |27| : the exact 
Ward-Takahashi relation for the low-frequency dynami- 
cal pseudo-spin susceptibility is obtained, at zero tem- 
perature, as 
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with 




FIG. 1: Schematics of our double-dot system: two dots are 
connected via tunneling t, and each dot is connected to a lead 
via V. 
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where p a dT (SJ T ) is the density of states (self energy) for 
an electron at the Fermi level in the dot r. 

We note that the second line of Eq. © is the well-kwon 
Korringa relation (2(| in the context of NMR relaxation 
theory, and the first line is its extension to the case hav- 
ing a finite energy-level splitting (2?]]. Since the static 
susceptibility can be calculated by the exact solution, we 
need to evaluate the density of states p a dT and the self 
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FIG. 2: Schematic description of the energy states in the 
double dots connected in series: (a) SU(4) symmetric case, 
(b) asymmetric case, (c) symmetric case in a magnetic field, 
(d) symmetric case in a strong magnetic field, (e) asymmetric 
case in a strong magnetic field. 



energy S^ r . Fortunately, this can be done by exploiting 
the Fricdcl sum rule. First recall that the phase shift 8° 
of an electron in the double dots at the Fermi level is ob- 
tained from the average number of electrons (n^ T ) m ^ ne 
double dots: 8% — ir(n dT ) (Friedel sum rule). Then, the 
density of states and the self energy at the Fermi level is 
given by 



ttA ' 
A cot 5° 



- dr 



(10) 

(ii) 



with the resonance width A due to the mixing V. Note 
that the electron number (n^ T ) can be evaluated by the 
exact solution. 

Combining all the above relations, we can compute 
the conductance at zero temperature: the quantities in 
the right hand side of Eq.lJSJ, p" dr , Yf dT , Xo P s(0), can be 
evaluated by means of the Bethe-ansatz solution of the 
SU(4) Anderson model H|. 



III. DOUBLE-DOT SYSTEM 

We study the conductance in several cases in our 
double-dot system, which are schematically shown in 
FIG. □ 



A. Charge fluctuations in symmetric double dots 

Let us start with the double-dot system shown in FIG. 
Uta), where the energy levels of two dots are same, which 
we refer to as the symmetric dots in this paper: there 
are 4 degenerate electron states including spin degrees of 
freedom. In this case, from the expressions 0, (JSJ and 
@ we write down the conductance in the absence of the 
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FIG. 3: (a) Conductance in the cases of the SU(4) and SU(2) 
symmetric double-dot systems as a function of the renormal- 
ized energy level e* . (b) conductance on log scale: we can 
see distinct exponential dependence between the SU(4) (zero 
field) and SU(2) (strong field) cases in the Kondo regime. 



magnetic field as, 



4tt^ 
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By computing static pseudo-spin susceptibility xj ps (0) 
by means of the Bethe-ansatz solution of the SU(4) An- 
derson model, we evaluate the conductance as a function 
of the effective energy level e* j2j| . The results are shown 
in FIG. |21 When the dot- level e* is above the Fermi 
level, the conductance is small, since the resonant tunnel- 
ing does not occur. As e* goes down through the Fermi 
level, the conductance is enhanced by the Kondo effect, 
which is analogous to an ordinary single dot case. How- 
ever, in contrast to the single dot system, for which the 
conductance is saturated in the Kondo limit with deep 
e* , it continues to increase exponentially. The increase is 
caused by the inter-dot charge fluctuations enhanced by 
the inter-dot "orbital" Kondo effect |22U23l . l24| . Since the 
static susceptibility Xo ps (0) is inversely proportional to 
the Kondo temperature Tk ~ exp(— A/e*), the conduc- 
tance has the exponential dependence like exp(— 2A/e*). 

Note that the ordinary spin Kondo effect and the inter- 
dot Kondo effect both emerge in the above SU(4) sym- 
metric case. Therefore, in order to see the above char- 
acteristic enhancement of the conductance more clearly, 
we consider an extreme case with strong magnetic fields, 
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where the spin Kondo effect is completeiy suppressed. 
Shown in FIG. E^b) is the conductance in strong mag- 
netic fields (corresponding to FIG. Eld)). We can see the 
enhancement of the conductance due to inter-dot Kondo 
effect with SU(2) symmetry. In this case, the correspond- 
ing Kondo temperature is given by Tk ~ exp(— 2A/e*), 
so that the increase of the conductance, ~ exp(4A/e*), 
is more significant in comparison with the zero field case. 
These results are indeed seen in log-scale plots given in 
FIG. |3 which clearly features the exponential depen- 
dence of the conductance in the Kondo regime. 



B. Symmetric double dots: magnetic- field control 

It is seen from FIG. [3] that in the Kondo regime with 
deep dot levels, the conductance in the SU(2) case (strong 
field) is larger than that in the SU(4) case (zero field). 
This implies that the conductance may be monotonically 
enhanced in the presence of a magnetic field. To clar- 
ify this point, we focus on the field-dependence of the 
conductance for the SU(4) symmetric double-dot system 
(shown in FIG. |2Jc)) in the Kondo regime. Following 
the way outlined above, we can derive the conductance 
in this case, 



G 
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(13) 



By exploiting the exact solution of the SU(4) Anderson 
model in the Kondo regime (so-called Coqblin-Schricffcr 
model) , we compute the conductance as a function of the 
Zeeman splitting Ez, which is shown in FIG.0J Also, the 
effective Kondo temperature Tk(Ez) is plotted as a func- 
tion of the Zeeman splitting on the log-log scale in FIG.E1 
Here we assume that the direction of the magnetic field 
is parallel to spin"f . It is seen that the magnetic fields en- 
hance the conductance, in contrast to the ordinary Kondo 
effect in a single-dot system. The inter-dot Kondo effect 
is caused by the degenerate energy levels in two dots, 
which still possess SU(2) symmetry in strong magnetic 
fields. Since the pseudo-spin susceptibility xl ps (^) in- 
creases with the increase of the field, thus resulting in 
the enhanced conductance. In strong fields, the effective 
Kondo temperature, which is defined by the inverse of 
Xo P s(0), is given by H|, 

T K (E z )/T K (0) ~ {E z /T K (Q))-\ (14) 

so that the conductance increases as 

G ~ (T K (E z )/T K (0))- 2 ~ (E z /T K (0)) 2 . (15) 

Here we note that the conductance for electrons with spin 
parallel (anti-parallel) to the magnetic field increases (de- 
creases). This effect might be utilized for spin-current 
control by using double-dot systems. The above results 
agree with those of Borda et al. obtained by the numer- 
ical renormalization group analysis [24|. 



It is to be noted here that the SU(4) Kondo resonance 
has been observed not only in a double-dot system j2^ 
but also in a single vertical quantum dot whose symmet- 
ric shape gives rise to SU(4) internal degrees of freedom 
[29j. Also, STM experiments on a Cr(OOl) surface have 
found the SU(4) Kondo resonance, where the degeneracy 
of d xz and d^, states gives the additional orbital degrees 
of freedom |30| . Our results are also consistent with these 
findings. 
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FIG. 4: Conductance as a function of the Zeeman splitting 
in the Kondo regime. We also show the contribution of the 
electrons with spins parallel (anti-parallel) to the magnetic 
field. Note that Tk = Tk(Ez = 0) is the Kondo temperature 
of the SU(4) Anderson model. 
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FIG. 5: The effective Kondo temperature as a function of the 
Zeeman splitting on log-log scale. 



C. Asymmetric double dots: gate-voltage control 

Next, we consider how the conductance is influenced by 
the energy-level difference between the two dots, which is 
controlled by changing the gate voltage of each dot. We 
study two typical cases in the Kondo regime: zero mag- 
netic field (FIG. [2b)) and strong magnetic fields (FIG. 
He)). 

From the expressions and (|S|l- lllfl . the conductance 
at zero field is written as 



G 



5E 2 



(16) 
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FIG. 6: Conductance as a function of the energy difference 
SE between two dots. We take T K = T K (E Z = 0) for the 
SU(4) double-dot case (zero field) and T K = T K (E Z = oo) 
for the SU(2) double-dot case (strong fields). 



We compute the conductance as a function of energy dif- 
ferences SE. We also study the conductance in strong 
magnetic fields, where the system is completely polar- 
ized, and the remaining inter-dot charge fluctuations are 
described by the SU(2) Kondo model subjected to the 
energy difference <5£^(23, ulj- The results obtained 
in both cases are shown in FIG. El In contrast to the 
magnetic-field dependence discussed above, the conduc- 
tance decreases monotonically as a function of SE. The 
decrease is caused by the suppression of the inter-dot 
Kondo effect in the presence of finite energy difference 
SE. 

We note here that the ordinary "spin" Kondo effect 
still persists even in finite SE, giving rise to the enhanced 
spin fluctuations. Although such enhancement in spin 
fluctuations may not be observed in transport proper- 
ties, it should show up if we observe the NMR relaxation 
rate in the double-dot system. For example, the in- 
dependence of the NMR relaxation rate is exactly given 
by the function shown in FIG. 0| for large SE, it is en- 
hanced as (SE) 2 . 



IV. TRIPLE-DOT SYSTEM 

We now discuss how the above method can be used to 
calculate the conductance for multiple-dot systems with 
more than two dots. Here, we deal with a triple-dot sys- 
tem, and then briefly outline how to extend the method 
to A-dot systems. We will see that the conductance ex- 
hibits some characteristic properties under the control of 
the gate voltage and the magnetic field. 

Let us consider a triple-dot system, for which three 
dots and three leads as arranged as shown in FIG. [7| 
Inter-dot tunneling t (intra- and inter-dot Coulomb re- 
pulsions) is assumed to be sufficiently small (large) here 
again. Therefore, one of the three dots can accommodate 
an electron thanks to strong intra- and inter-dot correla- 
tions. We focus on the Kondo regime, where the energy 
levels in the dots are sufficiently lower than the Fermi 




FIG. 7: Schematics of our triple-dot system: three dots are 
connected via small tunneling t, and each dot is connected to 
a lead via tunneling V. Inter- as well as intra-dot Coulomb 
repulsions are assumed to be sufficiently strong. 

level. We note that a similar but different triple-dot sys- 
tem has been proposed recently 31] and its symmetry 
properties have been discussed. 

In the second order in tunneling t, we can calculate the 
conductance between two leads in the triple-dot system 
by the exact solution of the SU(6) Anderson model, be- 
cause there are 6 available electron states including spin 
degrees of freedom in the three dots. In this case, we 
can still utilize the formula Eq.lJSJ to calculate the con- 
ductance, where the strong inter-dot correlations among 
three dots are incorporated via the inter-dot susceptibil- 
ity Xops between two dots through which electric currents 
flow. 



A. gate-voltage control 

Transport properties for the above three-dot systems 
depend on how the current is observed. To be specific, 
we change the gate voltage attached to the dot 3 with 
keeping the voltage in the dots 1 and 2 fixed, and 
observe the conductance between the leads 1 and 2 as 
well as between the leads 1 and 3. 

The computed conductance is shown in FIG. 02a) 
as a function of the energy difference SE between the 
energy level in the dot 3 and those in the dots 1 and 2. 
We set the sign of SE positive when the energy level in 
the dot 3 is higher than the others. 

Let us first observe the current between the leads 1 
and 2. It is seen that the conductance increases with the 
increase of SE (> 0). This increase is attributed to the 
enhancement of the inter-dot Kondo effect in the pres- 
ence of the energy deference, which is similar to that 
for the double dots in magnetic fields discussed in the 
previous section. At SE = 0, the current flows via an 
SU(6) Kondo resonance (i.e. 6-fold degenerate Kondo 
resonance). On the other hand, for large SE, the SU(4) 
Kondo effect is realized within 4 lower states in the dots 
1 and 2. This gives the enhancement of the inter-dot 
susceptibility between the dots 1 and 2, resulting in the 
increase of the conductance. Accordin g to the exact so- 
lution of the SU(6) Anderson model [32], the effective 
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FIG. 8: (a) Conductance as a function of the energy differ- 
ence. Here, is the Kondo temperature for SU(6) triple- 
dot systems, (b)the number of electrons for the asymmetric 
triple-dot as a function of the energy difference. The solid 
line is the number of electrons in the dot 1 or the dot 2 per 
spin and the dashed line is that in the dot 3. 



Kondo temperature for large SE is given by, 

T^ f) (SE)/T^ ~ (SE/T^y 1 / 2 , (17) 

and the corresponding conductance is 

G~6E/T$p, (18) 

where T^ 6) = T^ ff) (SE = 0) is the Kondo temperature 
for SU(6) triple-dot systems. 

On the other hand, for SE < 0, the SU(2) spin Kondo 
effect occurs in the lower 2 levels in the dot 3, while the 
spin Kondo effects in the dots 1 and 2 are suppressed 
because the number of electrons in the dots 1 and 2 de- 
creases (see FIG.[S{b)). Also, the inter-dot Kondo effect 
between the dots 1 and 2 is suppressed. As a result the 
conductance decreases when the current is observed be- 
tween leads 1 and 2. 

If the current is observed between the leads 1 and 3, 
distinct properties appear in the conductance. As seen 
from FIG.|HIa), for large \ SE\ (irrespective of its sign), the 
conductance decreases because the energy difference sup- 
presses the inter-dot Kondo effect between the dots 1 and 
2. Notice that around SE/T^ ~ — 1, the conductance 
has a maximum structure, where charge fluctuations be- 
tween the dots 1 and 3 are slightly enhanced. Any- 
way, the conductance exhibits behavior similar to that 
observed in the double-dot case under the gate-voltage 
control. 
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FIG. 9: Conductance as a function of the Zeeman splitting. 

B. magnetic-field control 

Let us now discuss how a magnetic field affects trans- 
port properties. For simplicity, we assume that the en- 
ergy levels of three dots are same (symmetric dots). The 
computed conductance between two leads under mag- 
netic fields is shown in FIG. The conductance in- 
creases as the Zeeman splitting Ez increases although 
the Kondo effect due to spin fluctuations are suppressed 
by the field. As discussed in the previous section, this en- 
hancement is caused by the inter-dot Kondo effect among 
three dots. For large magnetic fields, half of the internal 
degrees of freedom are quenched, so that the symmetry 
of the system changes from SU(6) to SU(3). As a re- 
sult, the SU(3) Kondo effect caused by inter-dot charge 
fluctuations is enhanced, and therefore the conductance 
is increased. The effective Kondo temperature in large 
fields is given as, 

T^ f) {E z )/Tf ~ {Ez/T®)-\ (19) 
and thus the conductance is enhanced like 

G~(E z /t£ ] ) 2 . (20) 



C. generalization to systems with more dots 

We can generalize our method to systems with more 
than three dots: a lead is attached to each dot, where 
the electrons feel strong intra- and inter-dot Coulomb 
repulsions. All the dots are connected to each other via 
small inter-dot coupling t. 

In similar manners mentioned above, we can calculate 
the conductance in such multiple-dot systems. The cal- 
culation can be done by using the formula Eq.J^J) in the 
second order in tunneling t, where all the correlation ef- 
fects are incorporated through the dynamical susceptibil- 
ity. We can use the exact solution of the SU(2iV) Ander- 
son model [25j for an TV-dot system. The conductance 
shows similar properties to those observed in the double 
and triple dots: if we change the gate voltage of the dot 
r, the conductance between the lead r and one of the 
other leads is generally suppressed, while it is enhanced 
otherwise. 
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In this paper we have assumed small inter-dot coupling 
t, and calculated the conductance up to t 2 . It should be 
mentioned that for a system with more than two dots, a 
Fano-type interference effect may emerge in higher order 
terms in t. This interference effect may give another in- 
teresting aspect of multiple-dot systems, which is to be 
studied in the future work. 



V. SUMMARY 

We have studied transport properties in the double- 
dot system connected in series that possesses not only 
intra- but also inter-dot Coulomb repulsions. It has been 
shown that the application of the Ward-Takahashi iden- 
tity enables us to use the exact solutions of the Anderson 
model for calculations of the conductance at zero temper- 
ature. We have clarified how the inter-dot Kondo effect 
affects the conductance under the gate-voltage control 
and the magnetic-field control. In particular, the con- 
ductance is decreased by the suppression of the inter-dot 
Kondo effect in the gate-voltage control, whereas it is in- 
creased by the enhanced Kondo effect in the presence of 
magnetic fields. The latter conclusion is consistent with 
the results of the numerical renormalization group. The 
method has also been applied to calculate the conduc- 
tance in multiple-dot systems including more than two 
dots. By taking a triple-dot system as an example, we 
have shown how the conductance is controlled by tuning 
the inter-dot Kondo effect. 

Naively, it seems not easy to observe the Kondo effect 
in multiple-dot systems (more than two dots) experimen- 
tally. We would like to mention, however, that the Kondo 
temperature in multiple-dot systems can be much higher 
than that in single-dot systems when the inter-dot repul- 
sion is relevant, as assumed in this paper. Therefore, if 



such multiple-dot system could be fabricated, the Kondo 
effect may be possibly observed even in multiple-dot sys- 
tems considered here. 

Finally a comment is in order on the ordinary "spin" 
Kondo effect in our multiple-dot system. We have 
focused on the inter-dot "orbital" Kondo effect in 
this paper, which directly affects transport properties. 
Concerning the spin Kondo effect, the impacts of the 
gate voltage and the magnetic field appear differently 
from the inter-dot Kondo effect, e.g. the magnetic 
field (gate-voltage difference) suppresses (enhances) 
the spin Kondo effect. If we use the dynamical spin 
susceptibility instead of the pseudo-spin susceptibility, 
the present analysis can be straightforwardly applied 
to low-frequency dynamics such as the NMR relaxation 
rate, which may be important to discuss an application 
to quantum-bits in quantum computation. In fact, the 
expression Eq.(JHJ) gives the NMR relaxation rate for 
the double-dot system, if the dynamical susceptibility 
is regarded as the spin susceptibility. It is of particular 
interest that the NMR relaxation rate in our multiple- 
dot systems can be controlled by the gate voltage, e.g. 
the difference in the energy-levels of double dots can 
enhance the relaxation rate. 



After the completion of this paper, we became aware 
of a recent preprint which deals with the SU(4) Kondo 
effect in a slightly different model|3j|. 
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